Hall coefficient in an interacting electron gas 
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The Hall conductivity in a weak homogeneous magnetic field, lu c t -C 1, is calculated. We have 
shown that to leading order in 1/cft the Hall coefficient Rh is not renormalized by the electron- 
electron interaction. Our result explains the experimentally observed stability of the Hall coefficient 
in a dilute electron gas in Si MOSFETs not too close to the metal-insulator transition. We avoid 
the currently used procedure that introduces an artificial spatial modulation of the magnetic field. 
The problem of the Hall effect is reformulated in a way such that the magnetic flux associated with 
the scattering process becomes the central element of the calculation. 

PACS numbers: 72.10.-d, 71.30.+h, 71.10.Ay 



I. INTRODUCTION 

In recent experiments on the two-dimensional dilute 
electron gases the renormalization of the Fermi-liquid 
parameters is found to be significant i - — These strong 
renormalizations are natural for a dilute gas with a large 
value of the ratio of the Coulomb energy to the Fermi 
energy. It was shown, however, that in a high mobility 
MOSFETs (metal-oxide-semiconductor field-effect tran- 
sistor) the density of electrons nn found from the Hall 
coefficient Rh = 1/nnec coincides with that obtained 
from the Shubnikov— de Haas measurements n to within 
a few percent [see Fig. 4 in Ref. |£| and the inset in Fig. 4 
in Ref. 0(b)]. This remarkable fact raises a question 
about the electron-electron (e-e) interaction renormal- 
izations of the Hall conductivity a xy in a small magnetic 
field, ui c t <C 1, where io c is a cyclotron frequency and 
t is the free path time. Notice that while the e-e in- 
teraction is very strong, in other aspects this system is 
rather simple, because the effect of the crystalline lattice 
can be safely studied in the effective-mass approximation. 
Apart from impurities, one may consider the electron liq- 
uid as if it is in a translational invariant background. 
Therefore the influence of the interaction on the Hall co- 
efficient should not be masked by unnecessary complica- 
tions. The issue of the Hall coefficient renormalization 
has not been discussed much in the literature despite its 
obvious importance. We have in mind the renormaliza- 
tion of the leading in T tr terms in the conductivity tensor, 
i.e., of the Hall conductivity a xy oc Br\ r , and the diag- 
onal term a xx = ne 2 Tt r /m*. Here r tr is the transport 
relaxation time, and B is a magnetic field; notice that 
Rh = B~ o xy j{a xx ) . 

In this paper we examine the effect of the e-e inter- 
action on the Hall coefficient within the microscopical 
theory. We avoid the currently used procedure that in- 
troduces an artificial spatial modulation of the magnetic 
field. Based on an idea of the flux of a loop in a dia- 
grammatic technique, we obtain a gauge invariant pro- 
cedure for the calculations in the presence of a homoge- 
neous magnetic field which makes it possible to conduct 
the analysis of a xy in a general manner. In particular, 
in the course of the calculations the translational invari- 



ance of the problem is maintained, and the usual rather 
involved procedure of extracting a constant B from the 
q — > limit of a singular expression for the vector poten- 
tial A(q) can be avoided. 

In this paper, we show that the Hall coefficient for elec- 
trons with the quadratic spectrum, e(p) = p 2 /2m, in the 
leading order in r tr , is not renormalized by the e-e inter- 
action, i.e., Rh = 1/nec. [In fact, the cancellation of the 
e-e renormalization in the Hall coefficient holds for any 
spherically symmetric dispersion e(|p|); see Appendix [B] 
for comments.] This result explains the stability of the 
Hall coefficient Rh to the e-e interaction demonstrated 
in Refs. and 0(b). Notice that the Altshuler-Aronov 
corrections to a xx and a xy teach us that there cannot be 
any general principle for the absence of renormalizations 
in Rh, because corrections to Rh of the order l/epTtr 
do exist £ Therefore the fact stated here about the can- 
cellation of e-e renormalizations in the leading term of 
Rh should be proved specifically. 

The renormalization of be also studied within 

the Fermi-liquid theory following the line given in Ref. 7 . 
The magnetic field B comes into a transport equation 
in a combination with the velocity of a charge carrier 
as (e[v p x TS]d/dp)5n p , where 5n p is the departure of 
the distribution function from equilibrium due to the 
applied electric field (see e.g., Ref. 8). In the pres- 
ence of the electron-electron interaction one may expect 
in addition terms like (e[v p x B]d/dp){f pp >5n p >), or 
(/pp'(e[ v p' x B ]d/dp')(fp>p"8n p »)), where (...) means 
the average over the Fermi surface. Here (f pp >Sn p >) is 
the response of the Fermi liquid to the departure of the 
distribution function from equilibrium. 7 This approach 
gives a xy renormalized by the e-e interaction which when 
combined with a xx yields the Hall coefficient Rh with no 
renormalization factors. It is, however, not clear that the 
derivation of the Hall conductivity in a weak magnetic 
field via the transport equation within the Fermi-liquid 
theory is accurate enough. For instance, consider the fol- 
lowing process: the electric field leads to the departure of 
the distribution function from equilibrium; the electron 
liquid responds to this Sn via the electron interaction; 
the magnetic field acts on the result of this response and 
transforms it into the Hall current involving again the e-e 
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interaction. The nonfactorized part of such two time in- 
teraction process (((e[v p / x"B]d/dp')f n p(p;p';p")5n p »)) 
corresponds to nonpair correlations, and is beyond the 
scope of the Fermi-liquid theory. (Some doubts about 
the derivation via the kinetic equation were expressed in 
the literature, e.g., Ref. 0, Sec. 7.12.) We will see, how- 
ever, that this contribution is small by the parameter 
l/f-FTtr- Therefore the present microscopical derivation 
justifies the applicability of the Fermi-liquid treatment of 
this problem. 

In Sec. II the technicalities of the calculation of a xy 
in a constant magnetic field are given. In Sec. Ill the 
Hall conductivity of the noninteracting electrons is con- 
sidered, and in particular a nontrivial problem of the Hall 
conductivity in a system with an arbitrary (nonshort 
range) impurity potential is studied. By studying this 
case we make the needed preparations for the analysis of 
the effects of the e-e interaction to the Hall conductivity 
[see Eqs. The problem of the Hall effect in 

a weak magnetic field is reformulated in a way such that 
the magnetic flux associated with the scattering process 
becomes the central element of the calculation. In Sec. 
IV the role of the e-e interaction is analyzed. There, the 
main contribution to a xy is given by Eq. I|35|) . The latter 
expression contains a product of two terms; one describes 
the scattering of the quasiparticles by impurities, while 
the other corresponds to the decoration of the current 
vertices by the e-e interaction. The derivatives over the 
momentum in Eq. I|35fl correspond to the vector product 
of the two coordinate differences, see Eq. (|23() . in the rep- 
resentation of the flux which leads to the skew action of 
the magnetic field. 

In the concluding section for completeness we discuss 
Rh using the Fermi-liquid theory^ In Appendix El we 
reproduce the known answer (Ref. ©) for the weak lo- 
calization corrections to o~ xy as an instructive example of 
the calculation within the new procedure. 

II. DIAMAGNETIC AND FLUX 
CONTRIBUTIONS IN A HOMOGENEOUS 
MAGNETIC FIELD 

To describe electrons moving in a magnetic field the 
magnetic vector potential A can be introduced through 
the extension of the momentum in the kinetic energy 
term: e(p) — » e(p — (e/c)A). Although the analysis 
presented in this paper is valid for any spectrum with 
a spherical symmetry (see Appendix [BJ) , the considera- 
tion in the main text will be limited to the case of the 
electron gas with a quadratic spectrum 

h = ^-( p --a) 2 , (1) 

2m V c / 

where m is the conduction band mass (i.e., m is not renor- 
malized yet by the electron-electron interaction) , and we 
use the fact that for low concentrations of the electrons 
only the quadratic term in e(p) is relevant. 



Notice that the homogeneous magnetic field demands 
special care as the Fourier components of the vector po- 
tential, A(q) = J d d re^ %<v A(r), are singular in this im- 
portant case. Let us take the vector potential in the 
Landau gauge, 

B = Bz; A x = 0, A y = Bx. (2) 

Then the Fourier components of A(r) are 

A y (d) = iB-^-8 d (q). (3) 

When inserted directly into diagrams the expression for 
A y (q) being highly singular leads to complications. In or- 
der to circumvent this difficulty it was proposed in Ref.lTfl 
to introduce a magnetic field modulated in space at some 
wave vector q in addition to a homogeneous electric field 
at finite frequency u>, and to study the (q,a>) response 
of the current. It has been argued in Ref. that in 
the calculation of o~ xy one should perform the limiting 
procedure as follows: first extract the q-linear contribu- 
tions related to the vector potential and arrange them 
into the gauge invariant combination B(q) = iq x A(q); 
next let q — > keeping B(q) to be finite; and, finally, 
perform the dc limit for the current response. The first 
step in this procedure requires combining different pieces 
into the gauge invariant combination, and this involves 
nontrivial cancellations amongst different diagrams. The 
outlined procedure is effective when a specific process is 
needed to be calculated (see, e.g., Refs. and0), but 
unfortunately it is unsuitable for a general analysis. 

The method proposed below treats the magnetic field 
at q = ab initio, and at no stage of the calculation 
does the vector potential appear explicitly. In addition 
to the apparent gauge invariance the developed technique 
shows clearly that although the Hamiltonian in the mag- 
netic field is not translational invariant, any measured 
quantity like conductivity can be calculated in an ex- 
plicitly transitional invariant manner. The calculation is 
based on the fact that the Green's function in a constant 
magnetic field can be represented in the formAi - — 

<7(ri, r 2 ,r) = exp [(ie/c)$(ri, r 2 )] 0(r x - r 2 , t), (4) 

where the phase $ is equal to 

$(r 1 ,r 2 ) = (r 1 -r 2 )A(^^±^) (5) 

and Q is the core Green's function which contains the 
information about the Landau-level quantization. The 
function Q is translational and rotational invariant and 
is also invariant with respect to different gauge represen- 
tations of the magnetic field B (see Appendix 1X1) . Notice 
that in the presence of e-e interaction the representa- 
tion (0J and (J5J is still valid. 

The rotational invariance of the core function Q is re- 
sponsible for the fact that the Hall current originates 
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either from the diamagnetic part of the current opera- 
tor, or with the participation of the phase terms $. No 
contribution to <r xy comes out from the core functions 
alone. This is because the calculation of a xy in the Kubo 
formula involves the averaging over the momenta of the 
nondiagonal current-current correlator (J x J y ). Since the 
functions Q are rotational invariant, the angular integra- 
tion of the type (p x . . .p') makes the core function's con- 
tribution to <T xy zero if the scattering by impurities is 
simple enough and does not itself produce the skew ef- 
fect on average. Then only the interference of the phase 
factor <E> or the diamagnetic current make it possible to 
achieve the skew action of the magnetic field to avoid 
the vanishing of the (J X J V ) current correlator after the 
angular integration. (Notice in this connection that the 
Lorentz force can be interpreted as a consequence of the 
phase accumulation in the magnetic field, Ref. ITS.) 

The diamagnetic part of the current operator, 
— (e 2 /mc)A, being nongauge invariant has to be com- 
pleted into a gauge invariant contribution in any calcu- 
lation of a measurable quantity. This completion comes 
out naturally here due to the differentiation of the phase 
factors in the Green's functions in the current operator, 

J v (ii) = (e/2rm) lim (V£ - V" )^ f (r^(r<) 

-(eVmc)^)^)^). (6) 

Consider the current vertex at the point connected to 
the rest of a diagram (e.g., for a current-current correla- 
tor) by a pair of Green's functions £?(ri, r$) and Q(r{, r 2 ). 
Then the diamagnetic part results in 

(eVmc^rO^n,!-,)^,^), 

whereas the first term yields 

(e/2m)[0(ri,r i )(-i)V r4 ^(r i ,r 2 )+iV r4 ^(r 1 ,r i )^(r i ,r 2 )]. 

To treat the phase factors in the Green's functions it 
is useful to split the A term into two half pieces and 
introduce the following relations: 

[-»V ri - (e/c)A(Ti)]g(T U r 2 ;T) = 

exp[(ie/c)$(ri,r 2 )] (7) 

x [-iV Ti + (e/2c) (r ( - r 2 ) x B] Q{n - r 2 ; r), 

and 

[»V ri - (e/cJAfafl^ri.rijT) = 
exp[(ie/c)$(ri,ri)] (8) 
x [-*V ri - (e/2c) (ri - r,) x B] g(n - r l5 r). 

The terms on the right-hand side that contain B explic- 
itly yield the diamagnetic contribution extended to gauge 
invariant combinations. Combining Eqs. Q and (JSJ) we 
obtain 

JtilZded = (ewe/4) [{Ti — r 2 ) x B — (n - r 4 ) x B 
x G(r 1 - ri)Q(ri - r 2 )R(r 1 - r/,r 2 - 17), 

(9) 



where R(ri — r_f,r 2 — Tf) is the remaining part of a cor- 
relator ending at a point 17. Notice that J^xtlnded ^ s 
also translation invariant. In Eq. @ the cyclotron fre- 
quency has been introduced in such way that its sign 
is that of the carriers: uj c = eB/mc. The phase factor 
$(ri,r 2 ) appearing in the Green's function ^(ri,r 2 ,r) 
may be rewritten as the integral over the straight line: 



$(ri,r 2 ) = (r 2 - ri)A 



ri + r 2 



A(r)dr. (10) 



Then collecting the phase factors from all Green's func- 
tions of a given loop of a diagram one will obtain the 
so-called "flux of a loop" of a diagram. [The phase fac- 
tors that appear in the right-hand side of Eqs. Q and 
IjSJ after passing through the differentiation operator join 
other such factors from the rest of the Green's functions.] 
It follows from Eq. i|10[) that the phase associated with 
each loop in a diagram is proportional to the flux, 



loop 



j> A(r)dr, 



(11) 



with the integral taken along the closed polygon corre- 
sponding to the loop. (The importance of the flux in the 
calculation of the Hall effect in the insulating state was 
emphasized in Refs . iTEl and ITil ) Notice that the polygon 
is oriented along the loop in the direction of the arrows 
of the Green's function. The magnetic flux through the 
polygon is gauge and translational invariant. It is equal 
to BS, where S is the oriented area of a polygon. De- 
pending on the situation one can decompose the polygon 
into the sum of different pieces. The decomposition can 
be done in a variety of ways, and a convenient choice 
may simplify the calculation. This is an advantage of this 
method. With the purpose of a general analysis of the 
Hall conductivity we will decompose the loop polygon in 
a diagram for the current-current correlator into the sum 
of triangles, and then use the fact that the flux through a 
triangle (ri,r 2 ,r 3 ) is equal to (l/2)B-[(r 2 -ri)x(r 3 -r 2 )]. 
The details of the decomposition procedure will be post- 
poned to the next section [see Eqs. (|20|l — (|23|l ]. where it 
will be shown how the flux contribution can be analyzed 
for the calculation of the Hall conductivity in a system 
with an arbitrary (nonshort range) impurity potential. 

Let us summarize. It follows from the symmetry of the 
problem that the Hall current may appear either from the 
diamagnetic part of the current operator, or through the 
phase factors of the Green's functions. Partially we have 
used the phase factors to extend the diamagnetic term in 
Eqs. 10— ©• The rest of the phase factors contributions 
can be organized in the form of the fluxes of the loops. 
The main consequence of this structure is that if one 
is interested in a xy linear in the external magnetic field 
only, a xy oc B, then it is enough for the calculations of 
a xy that the core Green's functions will be taken in the 
limit B — ► 0, when g coincide with the Green's functions 
in the absence of the magnetic field: £/b_>o = <?b=o- The 
latter fact enables us to reduce the calculation of a xy in 
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a small magnetic field to a diagrammatic problem that 
uses the Green's functions in the absence of the magnetic 
field. 



III. HALL CONDUCTIVITY FOR 
NONINTERACTING ELECTRONS IN THE 
PRESENCE OF A FINITE RANGE DISORDER 

We first illustrate the use of Eq. © by deriving the 
Drude Hall conductivity for the simplest case of short- 
range impurities. The conductivity will be calculated in 
the framework of the Kubo linear-response theory. The 
current-current correlator on the Matsubara frequencies 
is defined as 



U a/3 (iuj n ) 



e^ T (T T J a (T)J p (0)). (12) 



The analytic continuation from the discrete frequencies 
on the upper complex half plane to the real axes yields 
the retarded correlator H R Jui). This correlator is di- 
rectly related to the conductivity, 



o-a/aM = -n£a(g = 0,u;). 



(13) 



In the case of a short-range disorder the correlator Ijl2(l 
in the Drude approximation is equal to 



H xy (n,rf,r) = (e 2 /4m 2 ) lim {-iV x r , - {e / c) A x (v[)} + iV* - (e/c)A 



*(r,)) (-iV* 



(e/c)A»{T f ) 



- lim {i\7 r , - ( e /c)A»(i^)})g(r / -r < ,-r)exp[(*e/c)*(r / ,rO]g(r{-r / / ,r)ex P [(ie/c)$(4,r' f )] 



(14) 



Using Eqs. (JJJ and (JSJl, Eq. 114|l can be rewritten as 

U xy (n,r f ;T) = (e 2 /4m 2 )^lim_{-*V^ + (e/2c)[(rj - v' f ) x B] x } + *V* - (e/2c)[(r/ - r 4 ) x B] x ) ^-tV?, 
+ (e/2c)[( r/ - ri ) x B]„ + lim {zV« - (e/2c)[(rj - r^) x ByW, - r if -r)^ - r>, r). 

I 



(15) 



to momentum: 



Notice that when in Eq. i|15|l the gradients from the translational invariant one may use the Fourier transfor- 

J x vertex act on the coordinate dependent terms in the mation. Then the coordinate difference (r/ — r^) leads to 

J y part the corresponding contributions cancel out. The the differentiation of the Green's function with respect 
part linear in B arising from the J x vertex is equal to 

E^^r/jr) = ( e V/4m)[(r, - r/ ) x B] x 
x{HV« / g(r / -r i ,-r)]0(r i -r /) r) 
+ g(r f -r u -T)\iV« f g(r i -r f ,T)}}, 

and the other term from the J y vertex is 



(16) 



d d p 



(2^ (18) 

x[g A (p)d/d Py g R (p) - g R (p)d/dp y g A (p)] Py , 

and correspondingly 



U^( ri ,r f ;r) = (eV/4m) [(r f - r 4 ) x B] y 
x{[iV? t g(r / -ri,-r)]S(r i -r / ,r) 
+ g(r / -r i ,-T)[-iV*g(r i -r / ,r)]}. 



(17) 



(W (19) 

xp x [Q A (p)d/dp x g R ( P ) - g R {p)d/d Px g A {p)]. 



The summation over the frequencies and momentum in 
[Notice that Eqs. JEJ) and JEJ reproduce the struc- Eqs. JIHJ) and JTHJ) leads to the standard answer: a xy = 
ture of Eq. J^}.] Now the transition CJb^o = Ss=o can WcT<j xx . This part of the calculation is close to the one 
be performed in these expressions, and as all terms are presented in Ref. Il3i 
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Next we consider a nontrivial problem of the calcula- 
tion of the Drude Hall conductivity for noninteracting 
electrons in the presence of a finite range disorder, i.e., 
for nonshort-range impurities. To the leading order in 
1/epT the conductivity in this case is given by a set of 
ladder diagrams. To get a xy that is linear in the exter- 
nal magnetic field we consider separately the diamagnetic 
(extended) contribution, and the flux contribution that 
comes out from the phase factors accumulated by the 
Green's functions in a diagram. In the case of a finite 
range disorder the fcrmion loop does not degenerate to 
two retraced paths and the flux term becomes absolutely 
essential. (Naturally, we have to expand exp [(ie/c) ^i 00 p\ 
and keep the linear term only as we are interested in a xy 
that is linear in B.) 

The diamagnetic contribution is represented in Fig. ^ 
The cross means differentiation with respect to momen- 
tum. Only the Green's functions adjacent to the current 
vertex are differentiated. This contribution is obtained 
using the relations J7J), ©, and © similar to the case of 
the short-range impurities. 




b) 



FIG. 1: The diamagnetic contributions: (a) from the J x 
vertex, (b) from the J y vertex. 

To analyze the flux contribution to the Hall conductiv- 
ity we decompose the loop polygon in a diagram for the 
current-current correlator into the sum of triangles. The 
decomposition will be done as follows: one of the vertexes 
of a n-vertex polygon is chosen as a reference point (let 
us call it ri), and then ri is connected to all other ver- 
texes. In result the polygon is decomposed into oriented 
triangles all having the reference point as their vertex. 
(As an example consider a loop presented in Fig.|2) The 
flux through the polygon is therefore equal to the sum of 
the fluxes through each triangle: 

§loop = • [(r 2 - n) x (r 3 - r 2 ) 

+ (r 3 - r x ) x (r 4 -r 3 ) H ], 

which in turn is equal to 

$ioo P = • [(r 2 - ri) x (r 3 - r 2 ) + (r 3 - r 2 ) x (r 4 - r 3 ) 
+ (r 2 - n) x (r 4 - r 3 ) H ]. 

(21) 




FIG. 2: Flux contribution. 



Here we rewrote the vector differences (r\,- — ri) as [(vj — 

r-j-i) + (rj-_i ) + (r 2 - ri)]. Finally we add to $ /oop 

a zero flux expression that does not change its value, 

0=±B-(r n -r n _ 1+ r n _ a r 2 ^ 

+ r 2 - ri) x (n - r n ) . 

Now we are ready to state the general rule. After some 
reference point has been chosen the flux of the oriented 
polygon is given by 

lw 

where prime means that each pair or which 

is the same, enters the sum only once, and the pair is 
oriented with respect to the reference point in such a 
manner that its sign is "+" if we do not pass the reference 
point moving along the loop from r» to rj in the arrow 
direction. The sign is "— " otherwise. 

After the flux contribution has been decomposed in 
pairs each pair supplies the diagram with a vector prod- 
uct of the two coordinate differences. When Fourier 
transformed these coordinate factors lead to the differ- 
entiation of the corresponding Green's function with re- 
spect to the momentum components p x and p y . There 
are also situations when one or two of (r,-+i — ri) are ad- 
jacent to the current vertexes. Then the Fourier trans- 
formation leads, in addition to the differentiation of the 
vertex-attached Green's functions, also to the differenti- 
ation of the vertex momentum. 

Let the left-end current vertex be the p x vertex, and let 
the p y vertex be on the right end. The left-end located 
vertex will be chosen as the reference point with respect 
to which the orientation in Eq. I|23|l has been performed. 
First consider the terms represented in Figs.|31and0]when 
the vertices have not been differentiated. The difference 
between the two terms is that in Fig.[3]the differentiation 
d/dp x is closer to the left p x vertex, whereas in Fig. 0] it 
is d/dpy that is closer to the p x vertex. The ambiguous 
terms with d/dp x and d/dp y standing on equal distance, 
i.e., when they stand opposite to each other, are canceled 
out. Figure [S] represents schematically the flux contribu- 
tion involving the differentiation of the p y vertex at the 
right-end. In Figs.U2a) and (b) the d/dp x differentiation 
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FIG. 3: Flux contributions with no vertex differentiation 
involved and d/dp x is closer to p x vertex than d/dp y . 
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FIG. 4: Flux contributions with no vertex differentiation 
involved and d/dp y is closer to p x vertex than d/dp x . 
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FIG. 5: Flux contributions with p y vertex differentiation 
involved. 
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a) b) 
FIG. 6: Compensation of the false contributions of Fig. 



of the current operator: 



n 



(3,5) _ 
xy 



i (co c me 2 /2) [G A (p)d/d Px g R (p) 

g R (p)d/dp x g A (p)}j x d/dp y j y 



(24) 



(In the last expression we have omitted for brevity 
TJ2ie n /[ d V( 27r ) d ]> as wel1 as in E q s - <E3 and Jgj be- 
low.) In addition, Figs. Ufa) and (b) on one hand and 
Fig. 0Jc) and (d) on the other yield a similar combina- 
tion, 



= - i (cu c me 2 /2) \Q A (p)fl / d Vv Q R (p) 

- G^d/dpygAtpyjza/dpsJy. 



(25) 



acts on each place in the loop. It should be noticed, how- 
ever, that the momentum of the vertex originates from 
the gradients of the two Green's functions that are at- 
tached to the vertex [see, e.g., Eq. i|14[l]. Therefore the 
differentiations in Fig. [S] contain false nonflux contribu- 
tions when d/dp x and d/dp y act on the same Green's 
function attached to the vertex. The compensation of 
the false terms is given in Fig. where only the Green's 
functions closest to the p y vertex are differentiated. Fi- 
nally, in Fig. the contribution involving the p x vertex 
differentiation is represented. Here after cancellations 
amongst various terms combined into a full derivative 
only the two terms with the derivatives that act only on 
the Green's function closest to the p x vertex survive. The 
apparent asymmetry between Figs. |S] and compared to 
Fig. [7] is due to the fact that it was the p x vertex that 
was chosen here as a reference point for the orientation 
of the pairs, rather than the p y vertex. 

The total contribution. Let us collect all the terms to- 
gether. Notice the cancellations of Fig. with Fig. ^b), 
and Fig. □ with Fig. Then Figs. Ila),|3b), and^a) 
on one hand and Fig, ffic). E%d). and^b) on the other 
can be packed in combinations such that everything on 
the right from the p x derivative becomes a p y derivative 



Notice that finally we get a very simple rule: start with 
the derivative that stands at the left, and organize all 
the rest on the right as a derivative of the J y current; 
arrange the signs depending on whether the derivative 
is on the left-to-right or the right-to-left segments of the 
loop, and depending on whether it is a d/dp x or d/dp y 
derivative [the change of the sign is in accord with the 
orientation rules for pairs formulated for Eq. I|23fl ]. Re- 
markably, the diamagnetic and flux terms match each 
other to produce a rather simple form. (The Peierls sub- 
stitution for the hopping matrix elements in the presence 
of a magnetic field is useful to discuss the diamagnetic 
contribution on an equal footing with the flux term^ 
In the continuous limit both terms originate from the 
phase factor ~ie/c J rA in the path-integral formulation 
of the motion of the electron in the presence of the vector 
potential. 17 ) 

The Onsager relation, a xy (B) = a yx (— B), follows di- 
rectly from the rule formulated above. To get cr yx (— B) 
from a given contribution to a xy (B) one has to inter- 
change derivatives d/dp x and d/dp y and reverse B to 
compensate the change of the sign. 

Now we use a usual trick^fi and interchange in the 
IL x 4 y integration variables p x and p y . This leads to a 
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4 c 




b) 



FIG. 7: Flux contributions with p x vertex differentiation 
involved. 



convenient expression for H xy : 

n xy =iu c rJ- (g A ( P ) ^-g R (p) - g R (p)^-g A ( P 



dp x 

a d 

X 1 X dYy Jy ~ ^dVy^ 



dp x ' 



(26) 



The scattering by long-ranged impurities results in the 
renormalization of the current vertex by the ratio of the 
transport time r tr to the single-particle scattering time 



Jot 



Ttr{\P\)P a 



T 



m 



(27) 



Fortunately, in the combination (J x d/dp y J y — 
J y d/dp y J x ) the terms of the type d / dp y T tr {\p\) 
cancel each other out. Then the standard integrations 
in H xy lead us to a natural conclusion that the answer 
for the finite-range impurities differs from that one for 
the short-range case by the substitution r — * r tr : 



'xy 



U c Tt r Px 



(28) 



Correspondingly, it follows from this conclusion that the 
Hall coefficient is independent on the impurity range. 
[Here one point remains to be cleared up. In the presence 
of disorder the Green's functions acquire a self-energy 
part whose imaginary part is 1/t. Its real part shifts 
the chemical potential and hence is of no interest, but 
its dependence on the momentum near the Fermi energy 
may influence the result of the integration over momen- 
tum in Eq. (|26|l . However, in the absence of a special 
structure of the scatterers, like in liquid metals, the sen- 
sitivity to the energy of the electron state is very small 
when tFTtr 3> 1, and this effect can be ignored.] 

A comment may be in place here. It was important in 
the above derivation to keep the p derivatives in the form 
of combined expressions of the type d/dp y J y rather than 
to consider separate terms. Separately these terms are 
larger or even more singular than their sum. Let us see 
how it works in a somewhat pathological but instructive 
example of short-range impurities. In this case the dress- 
ing of the vector vertex by the impurity lines does not 
give any effect because of the momentum averaging. The 
p y derivative makes the vector averaging noneffective as 
it is illustrated in Fig. EI Each of the three diagrams pre- 
sented in Fig. [S] contains the singular propagator of the 






a) b) c) 

FIG. 8: Cancellation of the singular terms. 



diffusion ladder diagrams (diffuson) that does not vanish 
after the averaging. This does not lead to any complica- 
tions, however, as the sum of these singular terms van- 
ishes identically. Indeed, the expression corresponding to 
the ending block (i.e., separated by the last vertical im- 
purity line) when summed over the three diagrams turns 
out to be a derivative d/dp y {g R p y g A ) , where p is the 
momentum circulating inside this block. The integration 
over p forces the whole contribution to be zero. As a 
result only the diagram that does not contain vertical 
impurity lines survives in the case of the short-range im- 
purities as it is obvious when one considers these terms 
not separately but in the combined form d/dp y J y . The 
observation about the importance of keeping combined 
expressions in the calculation of a xy is of general charac- 
ter. Separate terms can be more involved than their total 
contribution, see, e.g., the discussion related to 9 2 E/<9£ 2 
following Eq. I|36[) . 

To conclude this section we notice that Eq. I|26[) repro- 
duces the structure of the iterative solution of the trans- 
port equation for a xy . Let us compare Eq. I|26[) with the 
term (e[v p x B]<9 / 'dp)5n p that appears in the transport 
equation as a result of the action of the Lorentz force and 
discribes the turn of the current. In the first bracket of 
Eq. (|26[) the derivative d/dp x yields the velocity v p , while 
d/dp y in the second bracket corresponds to the momen- 
tum derivative acting on the distribution function Sn p . 
The derivation in this section was limited to the leading 
order in \ jtpTtr- To get further, one has to study the 
dependence of f/rt,- on the flux. The leading term of 
that kind arises due to the Cooperon corrections and is 
discussed in Appendix IO 



IV. ELECTRON-ELECTRON INTERACTION 

We are interested in the current-current correlator, and 
therefore the basic process to be considered is a sequen- 
tial rescattering of electron- hole pairs. Diagrammatically 
this process is described by the particle-hole ladder sec- 
tions alternating with the irreducible amplitudes of the 
e-e interaction Y\ and T2- Amplitudes Ti and T2 differ in 
spin structure as it is schematically represented in Fig. [5] 
Examples of such amplitudes are given in Fig. I1UI 

A typical diagram for the Hall current correlator H xy 
includes the Xc-current vertex, several V amplitudes con- 
nected by the R ~ R, A — A, or R — A sections and, 
finally, the Jy-current vertex. In essence, the correlator 
Tl xy may be evaluated following the same line as in the 
case of the long-range impurity scattering. There is, how- 
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FIG. 9: The irreducible amplitudes of the e-e interaction I\ 
and T2. 





FIG. 10: Examples of the Vi and F2 amplitudes. 



and r; — > are used to form a new flux loop encircling 
the brick of the I\ amplitude in the direction opposite 
to the main flux loop. In Fig. Illf c) all the parts of the 
corresponding loop are given by the thick lines. 

To get the contribution to H xy from the main flux loop 
we shall follow the derivation given in the previous sec- 
tion. Consider first the situation when the left standing 
derivative acts on a Green's function in a particle-hole 
section rather than on the T amplitudes. Let this left 
standing derivative be the p x derivative. Then, the p y 
derivative is supposed to act on everything on the right 
of the p x derivative. The new element here is that the 
dashed lines added to restore the flux structure do not 
carry any Green's functions, but only the phase factors. 
Let us study the situation when the additional phase fac- 
tors arise from one of the Li amplitudes with the dashed 
segments r3 — > r 4 and rs — > Tq, as shown in Fig.^J Af- 
ter the decomposition procedure the corresponding flux 
contribution is equal to 



ever, a complication related to Ti insertions. In the case 
of noninteracting electrons we have dealt with a flux of a 
single loop extending from the J x to the J y vertex. An 
insertion of r 2 leaves the number of loops in a diagram to 
be the same, while the Ti amplitude splits the loop into 
pieces and breaks the single loop structure. We restore 
the construction with the main flux loop having in mind 
to use the basic results of the previous section, i.e., Eqs. 
PH l and P5|l. 

Consider an insertion of Ti that splits a loop at the 
points i, j, fc, I. To restore the construction with the main 
flux loop we attach the phase factor i(e/c)(rj — ri)A((rj + 
Yj)/2) to the segment — > Yj in Fig. Illf aL The anal- 
ogous operation will be done with the segment — » r/. 
Simultaneously the compensating phase factors with the 
opposite sign are attached to the lines going in the op- 
posite directions, Yj — > and 17 — > r^. In Fig. 1111 all 
these segments are denoted by the dashed lines with ar- 
rows. These dashed lines carry phases but do not repre- 
sent any Green's functions. Now we regroup the phase 





a) b) c) 

FIG. 11: The phase factors of the Fi amplitude. 



^dashed = ^ B • K r 2 - r l) X ( r 4 -r 3 ) , 2g . 

+ (r 2 - ri) x (r 6 - r 5 )]. 

Next, ^dashed can be rearranged as • [(r 2 — ri) x (r 4 — 
r$)+ (r 2 — ri) x (r@ — r 3 )], and for any structure of the 
Ti amplitude one obtains 

^dashed = • [(r 2 - ri) x (r 4 - r,- + r 5 ) 



+ (r 2 - ri) x (r 6 - r 4 + r. t r 3 



(30) 



It follows as a result of this rearrangement that after the 
Fourier transformation the operation of the p y differen- 
tiation will move all along the Green's functions in the 
segments r 3 — > r§ and r$ — > r 4 . 
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factors. It is possible because we are interested in a xy 
that is linear in B, and therefore we have to expand and 
keep only linear terms in the phase factors. The phases 
of the two dashed lines, — * Tj and — » r;, are used to 
obtain the flux of the main loop extending from the J x 
to the J y vertex [a fragment of this loop is presented in 
Fig. Illf b) where it is indicated by the thick lines] . The 
other two phase factors attached to the segments — > t.- l 



FIG. 12: The treatment of the dashed lines in the Fi am- 
plitude; the left standing derivative acts on a particle-hole 
section. 

Since Eq. I|30|) is valid for any Ti amplitude, we come 
to the conclusion that the dashed lines do not spoil the 
general structure of the discussed contributions, and ev- 
ery Green's function on the "vertical" lines, like r 3 — > Yq 
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and rs — > r 4 , is p y differentiated on the same footing as 
any other Green's function. Thus in the considered case 
the p y derivative indeed acts on everything on the right 
to the p x derivative, as it should be expected. 

Now we continue with the analysis of the case of the 
left standing p derivative that acts on one of the Green's 
functions in the particle-hole section. The sections R — R 
or A — A need not to be considered as each of them 
yields two expressions that cancel each other out. By 
the R — A section we assume a set of the impurity scat- 
tering ladder diagrams that modify a current vertex J 
by a factor T tr /r. In the limit when the external fre- 
quency u> — » only one R — A section should be kept in 
the current-current correlator H a p(u>), because each such 
section yields an additional power of oj. Differentiation of 
the Green's functions in the R — A section splits the set 
of the ladder diagrams into two pieces that modify both 
vertices J x and J y by r tr /r. As a result the contribution 
analogous to that represented in Figs. and reduces to 
Eq. I|24l) , and in the same way the contribution analogous 
to Fig. 21 reduces to Eq. I|25|) . Altogether we come back 
to Eq. <j2*rJ|l but with the current vertices and the Green's 
functions dressed by the e-e interaction. 

The terms that have been discussed untill now deter- 
mine the leading contribution to a xy cx t\ t . Its analysis 
will be completed in the end of this section, but first 
we turn to the other case when the left standing deriva- 
tive, let it be the p x derivative, is inside one of the in- 
teraction amplitudes T, while the p y derivative acts on 
everything to the right from this amplitude [an exam- 
ple is presented in Fig. EH a )]- As a candidate for a 
contribution of the order rf r to a xy in the discussed pro- 
cess consider the case when the p y derivative is applied 
to the R — A section standing on the right from T. This 
differentiation produces the singular terms (d/dp y Q R )Q A 
and Q R {d/dp y Q A ) that separately yield the contributions 
cx t\ t but in the leading order they cancel each other out. 
The latter fact is obvious because we may take the inte- 
gration over p y by parts and then the R— A section yields 
a contribution that is only cx r tr . A similar contribution 
appears when the R ~ A section is located on the left 
from the amplitude T that is marked by the differentia- 
tion. Arranging both the terms in a way that the R — A 
section becomes free from the differentiation the contri- 
bution to a xy in the discussed case can be presented as 



a xy cx LU c T tr J x (dq x f dp y - dq y f dp x )J y 



(31) 



Here the current vertices J XtV are dressed by the e-e in- 
teraction, but unlike J they do not contain the factor 
Ttr/r. The symbol r* [as well as the right-directed bar in 
Fig. I13f b)] indicates that the derivatives are ordered and 
the derivative over p acts on the right to the q derivative 
where q is an infinitesimal momentum in the current ver- 
tex J y as shown in Fig. ll4f a). The momentum q has been 
introduced to get opposite signs for the derivatives acting 
on oppositely directed lines in the current vertex. Since 
the parameters of the electron liquid change on the scale 




FIG. 13: A contribution to a xy that involves nonpair- 
correlation amplitudes. 



p + q/2 



p + q/2 




p-q/2 



p-q/2 




a) 



b) 



FIG. 14: The definition of the Fourier components of the 
current vertexes J and the interaction amplitude F. 



of the Fermi momentum this contribution is small com- 
pared to the leading one by the parameter l/ej?r tr . The 
reason of this smallness is that in the discussed process 
the skew action of the magnetic field may develop only 
on an electron wavelength, while for the leading term 
the skew effect develops on a free path length. Notice 
that the processes just discussed involve the derivatives 
of the r amplitudes that generate terms with the non- 
pair electron correlations f np that were mentioned in the 
Introduction. The discussed term may become signifi- 
cant if apart from the wavelength and the mean free path 
there is another length scale in the problem. For exam- 
ple, near the superconducting transition the fluctuations 
of the order parameter may contribute significantly and 
change the Hall coefficient. 

The last contribution to Tl xy that remains to be an- 
alyzed is related to the flux loops encircling the Ti 
amplitudes. Consider the loop encircling the ampli- 
tude ri(ri, r 2 , r 3 , r 4 ). The flux through the polygon 
ri,r 2 ,r 3 ,r 4 may be written as the sum of the fluxes 
through the two triangles: 



$n = B • [(n - r 2 ) x (r 4 - ri) 
+ (r 3 - r 4 ) x (r 2 - r 3 )]. 



(32) 



We define the Fourier representation of Fi(ri, r 2 , r 3 , r 4 ) 
as follows [see Fig. 114( b)] : 



(33) 



f d d p f d d p' f d d q 

x exp{i[(p + q/2)n - (p' + q/2)r 2 
+ (p' - q/2)r 3 - (p q/2)r 4 ]}r(p, p'; q). 

Comparing Eq. I|32|) with Eq. i|33|l we obtain that the 
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product Ti$ri when Fourier transformed yields the fol- 
lowing combination that enters the correlator H xy : 



ri* r 



B 



d 



i) 



i) 



<9q \ dp' dp 



ri(p, P ';q). (34) 



As the amplitude Ti is taken in the limit B = it is a 
scalar function of its momenta, and therefore the expres- 
sion 134|) vanishes. 

Summarizing the analysis presented in this section we 
conclude that the leading contribution to a xy is deter- 
mined by the following expression: 



n, 



m 

d 

OPx 



= ieB- 

2c \ T 



dp x 



(35) 



j —J 

d P y 



d P y 



The derivatives over the momentum in Eq. (|35() corre- 
spond to the vector product of the two coordinate differ- 
ences in the representation of the flux in Eq. I|23|) . The 
current vertexes J x ^ y are not sensitive to the renormal- 
izations by impurities due to 1/epT <C 1, and are deter- 
mined by the Ward identity*^ 



Ju 



m 



(36) 



where £(p, e) is a self-energy of the Green's function in 
the presence of the e-e interaction. Eqs. I|35l) and 1|36|) 
correspond to the solution of the transport equation for 
a xy for the interacting electrons [see the discussion of 
Eq. concluding Sec. III]. Due to the specific struc- 
ture of the last bracket in Eq. (|35|l the second derivative 
d 2 Yi/d£ 2 , which is not conventional in the Fermi- liquid 
theory, does not enter the final answer. Altogether we get 
the factor (l + <9E(p, 0)/<9£ p ) three times: two from the 
vertexes and one from the p x derivatives of the Green's 
functions in the square brackets. The integration over 
momentum in H xy can be performed keeping the Green's 
functions in the pole approximation: 

where r is the free path time of the quasiparticles that 
includes the renormalization by the e-e interaction, and 

m/m* = [1 + 0)/dt p ] I [1 - 9S(p F , e)/de] . (38) 

The latter relation is a standard one for the micro- 
scopical theory of the Fermi liquid^ All three fac- 
tors [1 — 8Yj{pf, e)/de]~ 1 created by the Green's func- 
tions in Eq. (|35|l together with the three factors 
[l + <9£(p, 0)/<9£ p ] can be combined into the cube of the 
physical combination m/m*. After integration over £ p 



one power of m/m* will be eliminated and ultimately we 
obtain 



(39) 



with uj* = eB/m*c and a xx — ne 2 T tr /m* . This result 
implies that no Hall coefficient renormalization devel- 
ops due to the e-e interaction in the leading terms over 
l/e F T tr . 



V. DISCUSSION 

We have shown microscopically that the Hall coeffi- 
cient in a weak magnetic field is not renormalized by the 
e-e interaction in the leading order in l/epT tr . Although 
the terms that are not conventional in the Fermi-liquid 
theory appear in the intermediate stages, they do not 
enter the final answer. This remarkable fact is a direct 
consequence of the skew structure which arises in the Hall 
effect because of the magnetic flux. 

The result about the cancellation of the Fermi-liquid 
renormalization corrections in the Hall coefficient differs 
from the one obtained previously in Refs. [lfj and [191 It 
is formidable to make a comparison with the analysis of 
these papers because the distinction of the impurity scat- 
tering amplitudes from the e-e interaction amplitudes 
was not performed explicitly. These amplitudes have 
different structure in exchange of the frequency, how- 
ever. Therefore any treatment lacking clear distinction 
between these amplitudes is potentially dangerous. The 
analysis of the Fermi liquid in Ref. [19| did not repro- 
duce the well established results of the transport theory, 
because of the confusing treatment of the frequency inte- 
grations induced by the e-e interaction amplitudes. On 
the other hand, in Ref. ^3 the- authors were mainly con- 
cerned with the impurity scattering. One can check, in 
fact, that their factor [1 + dE'(p, 0)/de(p)} 2 in R H dif- 
fers negligibly from unity in the case of noninteracting 
electrons scattered by a random potential. 

The absence of the renormalization in the Hall coef- 
ficient corresponds to the result that can be obtained 
within the phenomcnological theory of the Fermi liquid. 
For completeness we reproduce it here following Sec. 3.6 
in Ref. 0- In the Fermi-liquid theory the excited states 
of the interacting electrons are described by the gas of 
quasiparticles with the effective Hamiltonian written via 
the distribution function of the quasiparticles n p (r,i): 



(40) 



Here (p, r) are classical variables, and the interaction 
term / pp ' = J dr' /(pr, p'r') is assumed to be local. In 
the presence of the vector potential A(r,t) one should 
make a transition to the pair of the conjugate variables 
(P,r): 

e P (r) = e P _( e/c ) A (r); n P (r) = n P _ (e/c)A (r), (41) 
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where the vector potential A(r, t) represents both mag- 
netic and electric fields: B — V x A, E = —{l/c)dA/dt. 
In the phase space (P,r) the flow of the quasiparticle 
density rip{r) satisfies the Liouville equation: 

dn dn de dn de 

To make the dependence on the external fields explicit it 
is convenient to reexpress Eq. (|4*2|) in terms of n(p, r, t): 



the Altshuler-Aronov effect can be ignored, i.e., when 
(e 2 /ha xx ) ln(l/Tr) <C 1. Furthemore it follows from the 
analysis of Ref. |2£J that while as we show the leading 
term is robust, the temperature corrections to the Hall 
coefficient remain small up to the tempertures substan- 
tially lower than 1/r. Combined with this observation 
the present analysis gives an explanation of the stability 
of the Hall coefficient Rh to the e-e interaction observed 
m Refs. Hand 0(b) for the dilute electron gas in MOS- 
FETs not too close to the metal-insulator transition. 



dn ^ de dn de dn ^ dn 

dt dp a dr a dr a dp a " dp a 

e f de \ dn 
+ - — xB — =I coU . 
c \dp J a dp a 



(43) 



Next, it is useful to introduce the deviations from the 
global and local equilibrium defined as 

<5n p (r) = n p (r) - n°(e p ), 5n p (r) = n p (r) - n°(e p (r)), 



respectively. The two are related as follows 

dn° 



(44) 



5n p (r) = 5n p (r) - — £ / pp ,£ry (r). (45) 



With the use of Eqs. lO and GSl, Eq. lO yields after 
linearization 

dSn _ dn° d8n e , _,. dSn , Jn . 
- W +eE a — +v a ^r + - (vxB) Q — =/ co „ (46) 

where v a — p a /m* , and in the relaxation time approx- 
imation I co ii = —8n p /T tr - The current density is also 
expressed through the deviation from the local equilib- 
rium 



J a (r) = V* ^-5n p . 



(47) 



To study the Hall effect one should look for the station- 
ary, homogeneous distribution 8n p . The static limit of 
Eqs. H46f) and l|47|l are identical to the corresponding 
equations for noninteracting electrons, if m* — > m and 
8fi — > Sn. Therefore the transport coefficients are given 
by the free electron expressions with the substitution m 
to m*, and a xy = uj*T tr a xx . 

It is clear form this discussion that since the Fermi- 
liquid theory uses a local functional with the distribution 
function depending on the classical variables (p, r) it does 
not contain the nonpair-correlation contribution l|31|) . 
The latter describes the influence of the flux phase on 
the interaction amplitudes and, obviously, is beyond the 
scope of the Fermi-liquid theory. Fortunately, this term 
is small by the parameter l/ej?r tr . 

To summarize, we have proved that the Hall coeffi- 
cient is not renormilizcd by the e-e interaction in the 
leading order in l/e_FT. The result holds for not too 
low temperatures when the logarithmic corrections from 
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APPENDIX A: SYMMETRY PROPERTIES 

In this appendix the symmetry properties of the core 
Green's function G are presented. A useful discussion of 
the symmetry properties of the Green's function in the 
presence of e-e interaction can be found in the classical 
paper a 21 ii — which study the quantum oscillations in a 
magnetic field. 

1. Translation invariance 

A many-body system in a homogeneous magnetic field 
is described by the Hamiltonian 

^ = Ei(p-^)) 2 + I^mr l -r,|). (A!) 

In this case the vector potential depends linearly on the 
coordinates, and we introduce a matrix A^ v following the 
considerations of Ref. 



^(r) = A^r v . 



(A2) 



The conjugate vector potential is defined through its 
transpose A^(r) — A^ v r v in such a way that the gauge 
invariant relation A^ v — A^ v = Be^ holds. The genera- 
tors of the magnetic translations are defined as 



(A3) 



where the sum is taken over all electrons. Con- 
sider the operator of a finite translation S a — 
exp (— ia^T^). The Matsubara Green's function 
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^(ri,r 2 ;r) = — (T T ip(ri, r)^(r 2 )) remains unchanged 
under the canonical transformation ip —> ip = S~ 1 ipS a , 
as the generators (IA3(I commute with the Hamilto- 
nian lAlj l. Consider now this transformation explic- 
itly. Using the Baker-Hausdorff identity exp(A + B) = 
exp(A) exp(B) exp(— [A, B]/2) the S a may be decom- 
posed as 



S a =exp f-ia • P t 



x exp 



le 



t) exp ( 



i-a • A t 



j^p- _|_ 



(A4) 



where the Ptot is the total momentum and A to t is the 
sum of the conjugate vector potentials felt by all the par- 
ticles in the system. Then the transformed operator tp' 
becomes 



V>'(r) = -0(r — a) exp ( — aA(r) 



and correspondingly 

<,'t/Vi — 



~~ 

%j) T (r) = -0(r — a) exp ( aA(V 



(A5) 



(A6) 



Relations (|A5D and 1|A6|) imply that 

5(ri,r 2 ;T) = 0(>i - a,r 2 - a;r)exp ( — aA^ - r 2 

(A7) 

It follows then from the relation (|A7(I that the core 
Green's function defined in Eq. Q is translation invari- 
ant. In the presence of disorder the statement is still 
valid for averaged Green's functions. (The average over 
disorder generates translation invariant two-body "inter- 
action" terms when one introduces n — > replicas of the 
system, or uses other standard tricks.) 



2. Gauge invariance 

It has been shown in Ref. that Green's functions are 
transformed as 

exp (i-\x(ri) - x(r 2 )]) Gx{r\, t 2 ;t) = £/A+v x (n, r 2 ; t) 

(A8) 

under the gauge transformation A — > A' = A + V%. 
Clearly, the integral (ie/c) f* 1 Adr in the exponential 
prefactor in Eq. 0} changes exactly in a way as to cancel 
out the factor appearing in the left-hand side of Eq. I|A8() . 
This cancellation ensures the gauge invariance of the core 
Green's function. 



In this gauge the Hamiltonian (|A1|) takes explicitly the 
form which contains rotational invariant terms only. Sim- 
ilarly to S a the operators = exp (—i(j)L z ) can be in- 
troduced. As Sfj, commute with the Hamiltonian, the 
Green's function remains invariant under the action of 
S<p which rotate the coordinates ri and r 2 in £/(ri, r 2 , t). 
The phase factor (ie/c) f* 1 Adr is also unchanged in the 
course of rotation in the circular gauge, and therefore the 
core Green's function G should be rotational invariant as 
well. Since G is gauge invariant the statement holds for 
any gauge. 



APPENDIX B: SPHERICALLY SYMMETRIC 
SPECTRUM 

Transport properties in the case of a complicated band 
structure are nonuniversal. For this reason we are not go- 
ing to extend the present analysis of the renormalization 
of a xy by e-e interaction for an arbitrary spectrum. How- 
ever, the obtained result can be readily extended to Rh 
in the case of a general but spherically symmetric spec- 
trum e(|p|). (Such spectrum can naturally appear if the 
degrees of freedom in some energy shell are integrated out 
in the course of the Renormalization Group treatment.) 
In the presence of the magnetic field the corresponding 
part of the Hamiltonian is given by the gauge invariant 
extension of e(|p|) with a subsequent symmetrization re- 
quired to make the Hamiltonian Hermitian. The latter is 
needed since different components of the velocity opera- 
tor do not commute. To get a symmetrized extension of 
the kinetic term we expand e(|p|) in the Taylor series, 



IpI) 



/ J ^m.nPx Py • 
m,,n 



(Bl) 



The symmetrization procedure leads to 



e(|p|) =► £ e m , n (C« +n ) 1 ]T (p - -A 



PfH 



(p--a) ...( p -*a 

V c / P\2] \ c 



(B2) 



P[m+n] 



where the sum is taken over all incquivalcnt permutations 
P, in such a way that altogether there are m and n factors 
[p — (e/c)A] x and [p — (e/c)A] y respectively, and P[l] = 
x, y. The current operator is determined by the variation 
of the energy with respect to the vector potential: 



3. Rotational invariance 

We consider the transformation properties under rota- 
tions in the x — y plane in the circular gauge, 

1. 



-B x r. 



(A9) 



J M = -5(H)/5A l j i . 



(B3) 



Now we follow the main text after Eq. JHJ). Consider the 
diagrammatic element related to the current operator at 
the point after the contractions with the operators 
^(ri) and ij^ (r 2 ) have been performed. Using Eqs. I|B2|) 
and i|B3|) and integrating by parts we obtain for this el- 
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ement 



c / P(l-l) 



P I 

■ ■ • f-*V ri - -A) G(n - r.0 f-iV ri - -A) 

V c /P(i) V c /p(i+i 

•••(-£V 



£ A) G(r, - r 2 ). 

C / P(m+n) 



(B4) 



Here for each permutation all terms with the projection /i 
are varied to get the current ; correspondingly the sum 
^2' is taken over all I with P[l] = \x. In the absence of 
the vector potential the current operator yields de/dp as 
expected. To obtain the diamagnetic term in its extended 
form [see Eq. © above] it is needed to apply Eq. Q 
and to pass the terms [— iV — (e/c)A] through the phase 
factors using Eqs. Q and JSJ. The result is 



J™'" = e m ,n £ ^'(Crn+n)" 1 exp[(ie/ C ) (*(n , r<) 
p i 

+^,r 2 ))](-,V ri ^(r 1 -r i )xB) pM 

• ' ' f-*Vr, - |- (ri - r s: ) x b) G(ri - r<) 

V 2c ^P(i) 

X f-tV r< + 7T ( r i - r 2) X B ) 

V 2c /P(z+i) 



' • • -*V r( + — (r 4 - r 2 ) x B G{r t ~ r 2 ). 

V 2c / P(m+n) 

(B5) 

Here the gradients not acting on the Green's functions are 
canceled out due to the symmetrization procedure. At 
this stage the analysis of the main text can be extended 
up to the Eq. H35fl . where now 



J _ de(p) / [ 0S(p,6) 
The result still holds with 



(B6) 



(B7) 



It follows from Eq. (|B7fl that the renormalization correc- 
tions to Rh are canceled out. 



e'(p) fl + dZ(p,0)/dt; p \ eB 



P V 1 — 9^(pf, e)/dej c 



where <3> is given by Eq. J5J and C has the same symme- 
try properties as the corresponding function G in Eq. 
The latter fact can be seen from the Dyson equation for 
the cooperon in the magnetic field 

2-ktv 2lTTV J (C2) 

xg A (r u r')C(r',r f ) 

Using Eqs. (@J and (|C1|I one gets from Eq. I|C2|) an inte- 
gral equation for C, 



C(v u v f ) 



%-r/) , 1 



2tttv 



2-KTV 



£ A (r 4 -r)exp (2ie/c)$(r l; r , r f ) C(r ,r f ), 



(C3) 



where <J?(ri, r 2 , r 3 ) is the magnetic flux through the tri- 
angle with vertexes (ri,r 2 ,r 3 ): 



$(n,r 2 ,r 3 ) = ^B • [(n - r 3 ) X (r 2 - n)] . (C4) 



Equation (|C3|) contains invariant ingredients only, and 
this ensures that C has the same symmetry properties as 
G. One can readily derive the following equation for C: 



, „ 2eB 

D -iV x r 

2c 



C(r, w r , 



<5 d (r) 



27TT 2 ^ 

(C5) 

As it has been already explained in the main text for a xy 
linear in B it is enough to keep G and C in the limit 
B = 0. At the vanishing magnetic field C turns into the 
singular propagator 



1 



1 



2vrr 2 ^ taJ + DQ 2 



(C6) 



In this Appendix we will assume that the scattering 
potential is short ranged and that the electrons have an 
arbitrary but spherically symmetric spectrum e(p). Con- 
sider first the diagram presented in Fig. llSf a'). With the 



APPENDIX C: WEAK LOCALIZATION 
CORRECTIONS TO THE HALL COEFFICIENT 

The weak localization corrections originate from the 
quantum interference, and diagrammatically are known 
to be related to a set of the ladder diagrams in the 
particle-particle channel which is called "the cooperon." 
As a doubly charged object the cooperon in the magnetic 
field acquires the form 




"2 r 3 r 3 

a) b) c) 



FIG. 15: Weak localization corrections to u xy . The wavy line 
denotes the cooperon, the dashed wavy lines carry phases, the 
dashed line denotes the scattering by impurities. 



C B (ri,v f ) = exp [(»2e/c) $ (r^r/)] C(r 2 - r f ), (CI) 



use of Equation (0 the diamagnetic contribution from 
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the J x vertex can be written as 

™ / n 4 \ - e3B rrX- f dd P ( dd ® 

n xy (d = o, l « n ) = -s-rg J—J — 

dp x dp^ V dp v dVv j 

xC(Q)g R {- P + Q)g A {- P + qk i-p + Q) > 

(C7) 

where v a (p) — de/dp a is the velocity and d 2 ej dp a dpp 
is the inverse mass tensor. Together with the J y vertex 
contribution a symmetric combination 

[(d 2 e/dpl)(de/dp y ) 2 - 2(d 2 e/dp x d Py )(de/dp x )(de/dp y ) 
+(d 2 e/dp 2 y )(de/dp x ) 2 } 

arises. Exactly the same combination appears also in 
the calculation of a xy in the Drude's approximation. In 
the spherically symmetric case this combination reduces 
to (2/d)vpin~ 1 , where mT 1 — vp/pp. [The fact that the 
discussed combination contains the first derivative of e(p) 
only was important for the derivation of Rh = 1/ nec for 
the arbitrary spectrum e(|p|), see Eq. HB7L ] 

To evaluate the flux term, we split the phase of the 
cooperon into two equal parts indicated by the dashed 



lines in Fig. ll5f a). Then the flux of this diagram becomes 
a sum of the fluxes through the two triangles: 

$ a = *(r i ,n,r2) + $(r/,r 1 ,r 2 ). (C8) 

When the Fourier transformation is performed the coor- 
dinate differences in $ a lead to the differentiation with 
respect to momentum. To get a nonvanishing contribu- 
tion to a xy one of the differentiation should be applied to 
a current vertex. It turns out that the diamagnetic and 
flux terms are identical. Together they yield 

Sa a xy = 3cj c t5<t xx , (C9) 

where uj c = (eB/c)m , , and Sa xx is the weak localiza- 
tion correction to the longitudinal conductivity. 

For the terms presented in Figs. ll5f b) and (c) the dia- 
magnetic contributions are absent. The flux term of the 
diagram I15f b) may be represented most economically as 

$ 6 = $(r i ,ri,r 3 ) + $(r / ,r 2 ,ri) - 2$(ri, r 3 , r 2 ). (CIO) 

The advantage of this representation is that the first two 
fluxes are similar to the diamagnetic terms and do not 
contribute to a xy . The remaining flux term leads to 



n^(q = o,i« n ) 



.e 3 B ( 1 



V 2lXVT 



d d p f d d p' f d d Q t_\nAi snK, ^£ A (|-p + Q| 



(27r) d J (2n) d J (2tt) 

dg A ( P ') 



^v x (p)g A (p)g R (p)- 



d(- P + Q) x 



(Cll) 



x C(Q)g A (\-p' + Q|)^(|-p' + QD^^Vyi-p' + Q) 



The differentiated Green's functions are indicated by crosses in Figs. EP 3 ) an d (c). Similarly, the flux term of the 
diagram EI C ) is 



e 3 B ( 1 



2"KVT 



r d d p r d d pl r d d Q ggR 



R , ,^ R (|-P + Q|) 



-P + Q)* 



xC(Q)g A (\-pi + Q\)g R (\-pi + Q\) 



dg R ( P ') 
dp' Vy 



(C12) 



(-p' + Q). 



To proceed further, we rearrange the sum H xy = Ii b xy 
H xy to the following form: 

Ihrj^q = 0,iw n ) 



2-KVT 



dQ 



(27T) 



-MQ) 



%^ 2 [sUp)0A(p) + g 3 A (p)g R (p)} 



(C13) 



(2n) d d 

^L^Up)Ga(p) - Q\{p)Qr{ P )]. 



One should be cautious with the last integral as 



t{Q\q a - Q\Q R ) = 2it(g R g A f 



is an odd function in £. We have therefore to perform an 
expansion in £ to get a nonvanishing result. This leads 
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us to an integral 
2i f d(>irV 



-e {QrQat 



j d(yT 



£ 2 (S^a) 3 ■ 



(C14) 



Unlike the expressions discussed above, the integrand in 
Eq. I|C14|) depends on the specific form of the spec- 
trum, i.e., it is not universal. For e(p) = p 2 /2m and 
d — 2 the diagrams of Figs. lTBT b) and (c) yield Saxy — 



—u) c rba xx . Then the weak localization correction to the 
Hall coefficient, 



6R 



H 



Sa xy 



R 



H 



(C15) 



vanishes as it was first shown by H. Fukuyama^ 

The possibility of an economic representation of the 
flux is the essential advantage of the method. The expres- 
sions (|C7|) . IjClljl . and (|C12(1 have been obtained without 
any intermediate steps here. 
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